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Recently, a new model obtained from generalizing teleparallel gravity, named f(T) theory, is 
proposed to explain the present cosmic accelerating expansion with no need of dark energy. In 
this letter, we analyze the dynamical property of this theory. For a concrete power law model, we 



obtain that the dynamical system has a stable de Sitter phase along with an unstable radiation 
> " 



dominated phase and an unstable matter dominated one. We show that the Universe can evolve 
from a radiation dominated era to a matter dominated one, and finally enter an exponential 



expansion phase. 
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I. INTRODUCTION 

The models obtained from modifying general relativity theory are viable candidates 
to explain the current cosmic acceleration, which was first discovered from the supernova 
observations [lj and then further confirmed by many other cosmological tests, including 
the cosmic microwave background radiation [2( and the large scale structure [3j , and so on. 



[4] 



One such modification to general relativity is the f(R) theory (see |4| for recent review), 
where the Ricci scalar R in the Einstein-Hilbert action is generalized to an arbitrary 
function / of R. 

Recently, a new modified gravity to account for the cosmic accelerating expansion, 
named f(T) theory, is proposed by extending the action of teleparallel gravity [5|, |6] in 
analogy to the f(R) theory, where T is the torsion scalar. The teleparallel theory of 
gravity is built on teleparallel geometry, which was first introduced by Einstein to unify 
gravity and electromagnetism [5] and then was revived as a geometrical alternative to the 
Riemannian geometry of general relativity (7| . In teleparallel geometry, the Weitzenbock 
connection rather than the Levi-Civita connection is used. As a result, the spacetime has 
only torsion and thus is curvature- free. 

It has been demonstrated that the f(T) theory can not only explain the present cosmic 
acceleration with no need of dark energy 8|, but also provide an alternative to inflation 
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without an inflaton |9|, [10] • It therefore has attracted some attention recently. In this re- 



gard, Linder 11( proposed two new f(T) models to explain the present cosmic accelerating 
expansion and found that the f(T) theory can unify a number of interesting extensions 
of gravity beyond general relativity. We performed a statefinder diagnostic to these two 
models and also placed observational constraints on them from the latest data |12[. More 
recently, a reconstruction of f(T) theory from the background expansion history and the 



f(T) theory driven by scalar fields have been studied 13|, ll4|. In this letter, we plan to 
analyze the dynamical property of f(T) theory. 



II. THE f(T) THEORY 

In teleparallel gravity, the dynamical object is the vierbein ef , which has the property, 

e^ = S{, <#£ = %, (1) 

where e l is the inverse matrix of vierbein, i is an index running over 0,1,2,3 for the 
tangent space of the manifold, and (x, also running over 0, 1, 2, 3, is the coordinate index 
on the manifold. This vierbein relates with the metric through 

V = Vtj^ei , (2) 

where rjij = diag(—l, 1, 1, 1). 

As mentioned in the previous section, teleparallel gravity uses the curvatureless 
Weitzenbock connection, which is defined as 

rj, = ^e; = -e^e A . (3) 

From this Weitzenbock connection, one can introduce a non-null torsion tensor T a , 

T a = f A - f A . (4) 

[IV U/J, flV \ I 

The torsion scalar T in the action of teleparallel gravity is then given by 

T = VT% , (5) 

where 

S,T = \ {K^ + 5»T a » a - SZT*) , (6) 

and K^ u a is the contorsion tensor, 

K»\ = -\(T»\-T^-T a n. (7) 

For a flat homogeneous and isotropic Friedmann-Robertson- Walker universe described 
by the metric g^ = diag(—l, a 2 (t), a 2 (t), a 2 (t)) where a is the scale factor, one has, from 
Eq. ©, 

T = -QH 2 , (8) 



with H = aa~ l being the Hubble parameter. 

The action of f(T) theory is obtained by replacing T in the action of teleparallel gravity 
by T + f(T). Varying this action with respect to the vierbein, we obtain the field equation 
of f(T) gravity, which leads to the following modified Friedmann equation 

H 2 ~P-{- *&f T , (9) 

{H ] = 24/f/rr - 2 - 2/ r ' (10) 

where a prime denotes a derivative with respect to In a, subscript T, a derivative with 
respect to T, p is the energy density and P is the pressure. Here we assume that the 
energy components in the Universe are matter and radiation, thus 

1 . . 

P = Pm+pr, P = ~Pr , (H) 

where p m and p r represent the energy densities of matter and radiation, respectively. 

From Eqs. (jSJ [TO]) , we can define an effective dark energy, whose energy density and 
the equation of state can be expressed, respectively, as, 

Peff-^i-f^fr) (12) 

= f/T -f T + 2Tf TT + l^(f T + 2Tf TT ) 
Weff ~ (l + f T + 2Tf TT )(f/T-2f T ) ■ [ } 

III. DYNAMICAL ANALYSIS 

In order to analyze the dynamics of a general /(T) model, we rewrite the equations of 
motion as a dynamical system with the following dimensionless variables: 



where Q, r is the dimensionless energy density parameter of radiation. Using Eqs. OH [TO]) 
and the energy conservation equations of matter and radiation, one can obtain 

x' = -(2x + y) Z + 3 - 3X 9 - 3y , (15) 

2my — 2 + y 
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2/7 = 2my 2my-2 + y ' (16) 



,' = -Az - 2/ + 3 - 3x - % , (17) 

2my — 2 + y 



where a prime denotes a derivative with respect to In a and 



in 



h 



Defining r = — 2^- = -, one can express T as a function of y/x (or r). And then m 

J X 

can be expressed in terms of y/x. For example, the model, f(T) = a[(— T) p — f3] q , yields 
m(r) = (1 — q)r/2q +p — l. Thus, for a given form of f(T), the dynamical system given 
in Eqs. (TT5| [T6| H7|) becomes autonomous. 

Using x, y and z, Eq. (jUJ) and w e ff can be rewritten as 

^m = 3ij - 2 m =l-g-y-z, (19) 



x + y/2 — my 
(1 -y/2-my)(x + y) 



w eff = ~- 777^ Z^aTZ I ~ ' ( 20 ) 



where Q m is the dimensionless density parameter of matter. 

In order to analyze the dynamical properties of system given in Eqs. f|T5| [T6| ITT]) , we 
should firstly solve these equations with x' = 0, y' = and z' = 0. Here, besides two 
isolated critical points (denoted as Point A and Point B), we also get a continuous line of 
critical points, which is called as Line C: 

Point A: x c = 0, y c = 0, z c = 1 , (21) 

Point B : x c = 0, y c = 0, z c = , (22) 

Line C : x c = 1 — y c , z c = . (23) 

One can see that Line C is a straight line in the phase space. 
• Point A: radiation dominated point 



At this critical point, we have 

Q r = l , (24) 

which corresponds to a radiation dominated phase. Now we examine the stability of 
this point, which is determined by the eigenvalues of the linearized system. After some 
calculations, we find the eigenvalues at Point A, 



1, 2(1 - m ± Vl + Ira + m 2 - 2m') , (25) 

where m! = dm/ dr. This critical point is unstable because there is a positive eigenvalue. 

• Point B: matter dominated point 
Using Eq. (1T9|) . one has 

n m = i , (26) 

at this critical point. Thus, it represents a matted dominated era. Through the same 
calculation as that in Point A, we obtain that the eigenvalues of the linearized system at 
this point are, 

-1, 2 - 2m ± 2V1 + 2m + m 2 - 2m' . (27) 

For a successful cosmological scenario, this Point must be unstable so that the Universes 
can exit from the the matter dominated era. This means that the real part of one of the 
eigenvalues in the above expression must be positive. This may be possible for a given 
f(T) model as long as the model parameters satisfy certain conditions. 

• Line C: effective dark energy dominated era 
At this critical line, we have Q m = 0, Q r = and 

w eff = -1 . (28) 

Substituting Eq. ([9]) into f JTQj) and considering p = and p = 0, one can find easily that 

(H 2 )' = . (29) 

Thus, Line C corresponds to a de Sitter phase, if H ^ 0. The eigenvalues of the linearized 
system at this critical line are, 

-4, 0, -3. (30) 



It is easy to see that Line C is always stable. Thus, for a given f(T) model, the Universe 
finally enters a de Sitter phase. 

Now we consider a concrete power law model given in [8|, |ll| : 

f(T) = a(-T) n , (31) 



where a and n are two model parameters. In Refs. Ill, |l2j|, it has been pointed out 



that this model has the same background evolution equation as some phenomenological 
models 15Llla| and it reduces to the ACDM model when n — 0, and to the DGP model [l 



when n = 1/2. When n = 1, the Friedmann equation (Eq. (JOJ) ) can be rewritten as 
H 2 = 3 f^ a \ P, which is the same as that of a standard cold dark matter (SCDM) model 
if we rescale the Newton's constant as G — > G/(l — a). Thus, next, we will focus our 
attention on the case of n ^ 1. Let us note that, in order to be consistent with the present 



observational results, it is required that |n| <C 1 [8L llll. Il2|. 

Substituting Eq. ( |3Ti) into Eq. (jHJ), one can show that, when n ^ 1, the case p = 
(Line C) gives that the Hubble parameter is a non-zero constant, which corresponds to a 
de Sitter phase. Using Eq. (Tl8l) . we have m = — 1 + n and m' = 0. The eigenvalues at 
critical point B becomes 

-1, 2, -2n, (32) 

which means that critical point B is always unstable and its stability is independent of the 
value of model parameter n. Therefore, for a power law model, we find that, if n ^ 1, the 
Universe is able to evolve from a radiation dominated era to a matter dominated one, and 
finally enter an exponential expansion phase. In Fig. (CQ), we show the cosmic evolution 
with different initial conditions. It is easy to see that Xi and yi (f/T and —IJt at Oj, 
where a« is the initial value of scale factor) must be very very small to have a long enough 
period of radiation domination give the correct primordial nucleosynthesis and radiation- 
matter equality, and to ensure the appearance of a matter dominated phase, otherwise 
the Universe has unusual early behavior and evolves directly from radiation dominated 
phase to a de-Sitter one. Thus, for the power law f(T) model, the conditions for the 
Universe to evolve to a cosmic accelerating expansion and have usual early behavior are 



Xj >C 1, Hi -C 1 and n^l. Note, however, that in order to satisfy the current observations 



constraints, \n\ <C 1 is still required |S|, |ll|, Il2| 
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FIG. 1: The cosmic evolution for the case of power law model with n = 0.2. ctj is the initial 
value of the scalar factor. The dot-dashed, solid, and dashed lines correspond to the evolution- 
ary of curves of Q r , Vt m and the dimensionless density parameter of the effective dark energy 
respectively. In left panel the initial conditions are set as xi = yi = 10 -13 and zi = 0.98, while 
in right panel they are x^ = yi = 10~ 5 and Zi = 0.98. 

IV. CONCLUSION 



The f(T) theory, obtained from generalizing teleparallel gravity, is a new modified 
gravity capable of accounting for the present cosmic accelerating expansion with no need 
of dark energy. In this Letter, we analyze the dynamical behavior of the f(T) theory by 
assuming the existence of matter and radiation in our Universe. Two critical points (Point 
A and Point B), corresponding to a matter dominated phase and a radiation dominated 
one, respectively, and a critical line (Line C), corresponding to an effective dark energy 
dominated era, are found. We find that both Point A and Point B are unstable while 
Line C is always stable. Thus, the Universe can finally enter a de Sitter expansion 
phase, if H is nonzero at the critical Line C. For a power law model, the case n ^ 1 is 
considered since n — 1 corresponds to a SCDM model if we rescale the Newton's constant 
as G — > G/(l — a). The results show that, if n =£ 1, the final state of our Universe in 
the f(T) theory is an exponential expansion since H is a nonzero constant at Line C. In 
addition, we find that, to obtain the usual early universe behavior, it is required that Xi 
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and Hi (f/T and — 2/^ at Oj) should be very very small. Thus, for the power law model, 
the conditions to have a successful cosmological scenario are that Xj < 1, j/j C 1 and 
n/1. But, according to the results obtained in Refs. [8|, [ll|, ll2j, |n| <C 1 is required in 
order to be consistent with current observations. 
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